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We investigate in detail the effect of making the range of the "contact" interaction used in effective 
field theory (EFT) calculations of NN scattering finite. This is done in both an effective field theory 
with explicit pions, and one where the pions have been integrated out. In both cases we calculate NN 
scattering in the So channel using potentials which are second-order in the EFT expansion. The 
contact interactions present in the EFT Lagrangian are made finite by use of a square- well regulator. 
We find that there is an optimal radius for this regulator, at which second-order corrections to the 
EFT are identically zero; for radii near optimal these second-order corrections are small. The cutoff 
EFTs which result from this procedure appear to be valid for momenta up to about 100 MeV/c. We 
\Q ' also find that the radius of the square well cannot be reduced to zero if the theory is to reproduce both 

the experimental scattering length and effective range. Indeed, we show that, if the NN potential 
On ■ is the sum of a one-pion exchange piece and a short-range interaction, then the short-range piece 

must extend out beyond 1.1 fm, regardless of its particular form. 



O 



> 



ON 



x 



I. INTRODUCTION 



Effective field theory (EFT) techniques have been used successfully for many years to study problems in particle 
physics where a well-defined hierarchy of mass scales exists JlJ. In such problems one can, in principle, integrate out the 
short-range degrees of freedom, i. e. those corresponding to momenta above some separation scale, thereby obtaining 
a non-local effective action which will be applicable at energies well below this scale. The assumption is then made 
that all momenta in the problem of interest are small compared to the masses of the degrees of freedom which were 
integrated out. It follows that a momentum expansion of the action may be made. This leads to calculations which 
are organized in terms of the power of momentum which a given diagram contributes — the so-called EFT expansion. 

This is possible because power-counting arguments indicate that a given loop diagram only contributes at a definite 
order in the EFT expansion. Naturally, these loop diagrams are divergent, and must be regulated if they are to yield 
finite results. These results must then be renormalized. The theory as a whole is not renormalizable, but at a given 
order in the EFT expansion there are only a finite number of possible counterterms and so the theory does have 
predictive power. Hence, once the order of an EFT calculation is fixed, the predictions of that theory are independent 
' of the regulation scheme chosen. 

Several years ago, Weinberg suggested that the techniques of EFT could be modified to study low-energy problems 
in nuclear physics, including nucleon-nucleon scattering ^|,^). He pointed out that the fundamental difference in the 
NN scattering case is that the EFT expansion must be applied to the potential — not the amplitude, as in most other 
problems. If this is done then the potential which is to be used in the Schrodinger equation contains a delta function 
interaction. Direct calculation with such an equation is impossible. Regulation and renormalization of the results 
must be performed before physical predictions can be extracted. The question is how this is to be done. 

So far two approaches exist in the literature. The first involves formally iterating the divergent interaction, and 
then renormalizing the resultant amplitude. This was the approach used by Weinberg More recently, it has been 
applied by Kaplan et ai, who iterated the contact interaction formally and then used dimensional regularization and 
the MS renormalization scheme to calculate NN scattering in the 1 So channel [Q . The second approach is to introduce 
a regulator into the equation before iterating the potential to generate the t-matrix. This was the method adopted 
by van Kolck et ai, who cut off all integrals in momentum space in studies of NN scattering |5|. In such an approach 
renormalization is performed by adjusting the coefficients of the potential to fit the NN scattering data. This can 
be done for several different values of the regulator parameter. The sensitivity of unfit physical observables to the 
regulator parameter may be regarded as a measure of the validity of the regulation. If an EFT solution to the NN 
problem can be constructed by one of these two, or some other, means, it may then be used as input in nuclear physics 
applications of EFT, e.g. to pion-deuteron scattering the three-nucleon system [Q, the pp — > ppir process ^,|), 
and pion photoproduction on the deuteron Jlp| l . 

One might think that in order to establish the connection with the original field theory the regulator in the second 
approach must ultimately be removed from the problem, i.e., that the range of the regulated contact interaction must 
eventually be taken to zero. However, it has recently been pointed out that the use of zero-range interactions which 
are defined as the limits of short-range ones as the range is taken to zero leads to peculiar results. In particular it 
has been shown that EFTs either with or without pions cannot contain a repulsive zero-range interaction (defined in 
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this sense) which has any effect on the observables in the renormalized theory ]ll]J12| . Thus all non-trivial zero-range 
interactions must be attractive. Furthermore, in theories without explicit pions (i.e. where the pions are integrated 
out) it was shown that if the range of the interaction is taken to zero then the phase shifts obey: 

-±(kcot(5(k)))<0. (1.1) 

A corollary is that potentials of sufficiently small radii cannot reproduce the NN 1 So scattering length and effective 
range. In this paper we show that similar constraints apply when pionic degrees of freedom are explicitly included. 

While these results do not invalidate the conclusions of Ref. they do show that the approach used there is not 
equivalent to solving a Schrodinger equation with an interaction of range R, and renormalizing the coefficients as R is 
adjusted. In other words, for this problem, regulating the potential and then iterating it is not equivalent to iterating 
it formally and then using dimensional regularization and MS renormalization. 

This leads us to investigate in detail what happens if the range of the "contact" interaction is made finite before 
iteration, as was done in practice in Ref. 0. (The use of similar "cutoff" field theories has been advocated for 
some time by, for example, Lepage While such an approach apparently violates some of the assumptions of 

Weinberg's power-counting argumentsp] , one trusts that for a wide range of regulator parameter values the results of 
the renormalized theory are not particularly sensitive to the exact value of that parameter. Indeed, a virtue may be 
made of the necessity of keeping this parameter finite. The freedom to choose it may be exploited so as to minimize 
the error in truncating the EFT expansion. That is, a range may be chosen for our short-range interaction which 
results in the renormalized coefficient of the second-order term in the EFT expansion being zero. For other values of 
the short-range interaction range the value of this renormalized second-order coefficient will, of course, be non-zero. 
The viability of such an approach must be judged by testing how sensitive observables are to the range chosen for the 
regulated contact interaction. 

In this paper we explore this issue for two EFT- motivated treatments of nucleon-nucleon scattering in the Sq 
channel In Section |l| we investigate a Lagrangian with only contact interactions (i.e. where the pions have 



been integrated out), while in Section III we look at an EFT with contact interactions and explicit pions. In both 
cases we never use a true contact interaction but instead always solve the Schrodinger equation keeping the regulator 
parameter for the contact interaction, R, finite. The particular form of the regulator chosen here is a square well. In 
both EFT treatments we begin by calculating at zeroth order in the EFT expansion. At this order, the potential in 
the Schrodinger equation is either a square well alone, or the sum of a square well and one-pion exchange. For any 
given radius R the strength of the well can be adjusted to match the 1 So NN scattering length. However, by exploiting 
the freedom to choose R, we can minimize the second-order corrections to the EFT by choosing a radius which results 
in the total zeroth-order EFT potential reproducing both the scattering length and the effective range. The resulting 
value of the well radius, R, is the "optimal" one, since it minimizes the second-order EFT corrections to the potential. 
For an arbitrary regulator radius R, going to second order in the EFT means adding terms to the potential which 
correspond to derivatives of our "contact" interaction. The coefficients of these new terms are fixed by demanding 
that at different well radii the scattering length and effective range are still reproduced. At the optimal radius these 
coefficients are, by definition, zero. For radii close to optimal, it is found that the coefficients can be renormalized 
successfully. However, we will show that there is a lower bound on the radius for which this renormalization can be 
done. Below a certain value of R it is impossible to adjust the coefficients to reproduce low-energy 1 So NN scattering 
data. This explicitly demonstrates, for the specific case of a square-well regulator, that the theorem of Ref. 
prevents one from taking the limit R — ► 0. More generally, the result of Ref. |Q leads to an absolute lower bound 
on the radius R that may be used if the effective field theory is to correctly predict the effective range. This lower 
bound is completely independent of the particular regulator chosen. In the EFT without explicit pions it is a simple 



matter to show that the absolute lower bound is R = 1.3 fm. In Section IV we use a modified version of the argument 
presented in Ref. fll2| to show that the smallest possible i?, even when explicit one-pion exchange is included, is larger 
than 1.1 fm. 

The sensitivity of the results to the range of the short-range potential is explored in two ways. First, we examine 
the phase shifts for different values of R. Agreement is expected at low momenta since we have fit the scattering 
length and effective range. Thus the phase shifts should be identical until the point where the fourth-order terms in 
the effective range expansion become significant. The question is at what momentum these fourth (and higher) order 
terms become important for different well radii. Secondly, we examine the magnitude of the second-order terms in 
the t-matrix for potentials of various radii. The hope is that for radii close to optimal the second-order corrections 
are small for moderate momenta. This would indicate that if short-range potentials of these radii are used then the 
EFT expansion stays under control. 
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II. AN EFFECTIVE FIELD THEORY WITHOUT PIONS: THE SQUARE WELL AND DERIVATIVES 



As observed in Refs. if one integrates out all the exchanged mesons in the NN interaction then an EFT is 

obtained which consists solely of contact interactions and derivatives thereof. In the 1 5o channel the only pieces of 
this Lagrangian which contribute (to second order in the EFT expansion) are: 

C = N f id t N - N^^—N ~ ~C s (iv"+jV) 2 - ]-C t {N^cjN) 2 - jC 2 (N^V 2 N){N^N) - ^C 2 ((V 2 N^)N)(N ji N). (2.1) 

If the interaction terms here are used, as advocated by Weinberg, to generate a potential for use in the Schrodinger 
equation, then at second order in the EFT expansion the result is a Schrodinger equation which we will write with a 
general non-local potential, 
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V>(x) + / T/ fl (x,xXx')dV = £ty( x ), 



(2.2) 



where the subscript R indicates that V must be regulated in order to make the Schrodinger equation meaningful. The 
regulated non-local potential is: 



Vk(x,x') = (Co + C 2 (~V 2 ~ V' 2 ))^(x - x')^ 3) (x). 



(2.3) 
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Here Co = C? — 3Ct, and S R is a function which tends to a delta function as R — > 0. In this work we use a square 

(3) . 

3e(ij-|x|) 



(3) 

well of width R for the function 5 



4 3) (x) 



4nR 3 



In the I = partial wave, the radial equation obtained from the Schrodinger equation (2.2) is then 
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(2.4) 



(2.5) 



The solu tions to this equation in the two regions r < R and r > R must be matched at r = R. By integrating 
Eq. (2.5), we see that the boundary condition at r = R involves a discontinuity of the derivative of the wavefunction, 



due to the derivative terms in the potential: 



3C-AA dU 



4ttR 3 J di 



dU 
dr 



(2.6) 



where dU I dr \± d enote the one-sided derivatives at r = R. 

From Eq. (pit) an expression for the phase-shift 5 at momentum k = {ME) 1 / 2 may be derived. It is: 



k cot 5 — 



BKCot(nR) + kta,n(kR) 
1 - BfcoUnR) tan(fci?) : 



where 



B = 1 



3C 2 Af 
AttR 3 ' 



h 2 _ 3C M 
i , 3C 2 M 



1/2 



(2.7) 



(2.8) 



Expanding Eq. (2.7) in powers of k 2 , we find that the scattering length, a, and effective range, r e , of this potential 
are given by 



a = R 



R 



Bx cot x 



(2.9) 



R(R — a) 



2irR i (1 — xtanx — x cot x) 
(2ttR 3 + 6C 2 fi)x 2 a 



(2.10) 
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where 



R 



-3C M 



4vri? 3 + 6C 2 M 



1/2 



(2.11) 



As we can see from Eqs. (2.9) and ( 2.10) ) above, a; is a more natural variable for this problem than Co- Consequently, 
from now on we will generally use x and C2 as variables. 

The values adopted for the NN 1 Sq scattering length and effective range here are: 



-24 fm; 



2.7 fm. 



(2.12) 



We will examine how x and C2 must be ren ormali zcd a s R is varied in order to reproduce these values. We will 
begin by setting C2 = and solving Eqs. (2.£) and ( 2.10 ) numerically for R and x. We will refer to the result as the 
"optimal" well since it is the well which minimizes the second-order corrections in the EFT expansion. We will then 
vary R from the optimal value and solve for x and C2 to determine over what range of R we can successfully obtain 
the actual scattering length and effective range using this potential. Finally, we will examine the phase shifts and the 
second-order corrections to the t-matrices for these renormalized potentials and compare our results with the phase 
shifts extracted from the experimental data. 

Before we begin this procedure note that Phillips and Cohen (l2) have recently shown that a bound for dS/dk 
derived by Wigner [|15| yields the following constraint on r e : 



R 3 



d2 

-I « . 

a 6a z 



(2.13) 



Eq. (2.13) provides an absolute lower bound on the size of a potential which can reproduce a given scattering length 
and effective range. In the 1 S'o channel for NN scattering, this lower bound is 1.3 fm. In fact, for this potential, we 
will find that the lowest value of R for which Eqs. (2.9) and (|2.10| ) can be s olved is considerably larger than this. 

To determine the optimal width of the square well, Eqs. ( |2.9D and ( [2.10] ) are solved with C2 = 0, with the result 
that the optimal width is 2.6 fm. The corresponding value of x is 1.5. These values of x, C2, and R are equivalent to 
Co = —5.1 fm 2 . A plot of the phase shift as a function of k for these values of R, x and C2 is shown in Figure |. We 
confirm that, as we expect from Levinson's Theorem for a potential with no bound states, the phase shifts at zero 
and very large energy are equal. 
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FI G. 1. A plot of phase shifts in radians versus momentum for the optimal well. The phase shifts are calculated from 
Eq. (§/7|) with R = 2.6 fm, x = 1.5 and C 2 = 0. 



We now proceed to examine the behavior for other values of R. When Eqs. (2J)) and ( [2.10| ) are solved numerically 
for x and C2 as R is varied from the optimal value, we find the behavior shown in Figures |2j and [| As the radius is 
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reduced, we find that x tends to ir/2 and C2 diverges as we approach a radius near 1.7 fm. The fact that x is tending 
to a finite quantity while C2 diverges suggests that Cq is also diverging. For large R, C2 becomes large and negative. 
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FIG. 2. A plot of the renormalized values of x, defined as in Eq. (2.11), versus the radius, R, of a square well potential 
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FIG. 3. A plot of the renormalized values of the second-order coefficient in the EFT expansion, C2, versus the radius, R, of 
a square well potential. The value of C2 diverges near R = 1.7 fm. 



Armed with renormalized values of x and C2 for several different radii, we can use Eq. ( p.7[ ) to examine how the 
phase shifts are affected by changing R. Figure || compares the phase shifts for several values of R with those obtained 
at the optimal value of R. Since we have fit the scattering length and effective range, the phase shifts will agree for all 
radii until the fourth order terms in the effective range expansion become important. While this occurs at quite low 
momentum for the 5 fm well, when the radii are close to optimal the agreement persists to a much higher momentum, 
k « 0.5 fni . This is roughly what we expect since the radii of the relevant wells are themselves around 2 fm. Since 
this is a theory without explicit pions this radius can be thought of as being intrinsically that of one-pion exchange. 
Thus one might hope that when pions are explicitly included in the theory the radius of the short-range potential will 
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decrease significantly. 

For non-optimal wells the phase shifts always go to positive or negative infinity as the momentum gets large. 
However, this does not contradict Levinson's Theorem since these potentials are non-local. 
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FIG. 4. A plot of the phase shifts in radians versus momentum calculated from Eq. (2.7) for square wells with different radii. 
The phase shifts for the optimal well, shown in more detail in Figure [l], are shown here in bold for comparison. 

Next we compare these phase shifts with the experimentally determined values. Figure |5| compares the data from 
the Nijmegen partial wave analysis Jll| with the phase shifts from the optimal well and the phase shifts from the 
effective range expansion. Unfortunately, we find rather poor agreement with the data for our optimal square well. 
Contrastingly, the effective range expansion matches the data surprisingly well, indicating that the shape parameter 
for NN scattering in this channel is quite small. Of course, we could adjust all three parameters R, Co, and C2 in 
order to improve the agreement of our result with the data, but this would not be in the spirit of the calculations we 
are pursuing here. While the result displayed in Fig. |5| is not particularly promising for the approach espoused in this 
paper, we will see that the explicit inclusion of pion exchange will lead to a considerable improvement in the fit to 
the experimental data. 
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FIG. 5. A comparison of calculated and experimentally measured phase shifts. The solid line in this plot is the phase shift 
as a function of momentum for an optimal square well. The dashed line is the effective range expansion to second order. The 
dots are data points from the Nijmegen np partial wave analysis in the 1 5'o channel [ fl6| . 

Having examined the behavior of this system to second order in the EFT expansion, we now wish to examine just 
how significant the second order, i.e. derivative, terms are. For a given value of R, we may examine the fractional 
difference between the on-shell t-matrices at zeroth and second order in the EFT expansion, 

AT* fc ee R 1 > R { ' . 2.14 

Tf (fc) 

(2) 

Here, the second-order t-matrix, T# has already been calculated: 



where fccot<5 is given in Eq. (2/7). On the other hand, is the zeroth-order EFT amplitude for a square well of 
radius R. Hence, at a given R we h ave Ci = and the only free parameter is Cq. Consequently, we must renormalize 
differently. In this case we use Eq. fl2.S|) to fix x by demanding that we reproduce the x Sq NN scattering length. This 
gives the renormalization condition 

R 

XqCoIxq = . (2-16) 

R — a 



is then found by using the values x — xq and C-2 = in Eq. (2/7) for fccot<5. 
For the optimal width well, AT is zero, by construction. For non-optimal widths, Figure ^| shows AT^(fc) for several 
values of R. We note that, as we would expect, in the zero-energy limit the difference is always zero. The rate at 
which AT* grows with energy is dependent on how much the radius of the well differs from that of the optimal well. 




FIG. 6. A p lot of the fractional difference between the zeroth- and second-order t-matrices versus momentum, calculated 



from Eq. (2.14), for several values of R, the radius of the square well. For the optimal well with R = 2.6 fm, this quantity is 



identically zero. 

To summarize, we have shown how the coefficients in the second-order EFT expansion must be renormalized if we 
wish to reproduce the 1 S'o NN scattering length and effective range with a square well potential of a given radius. In 
particular, we have demonstrated that there is a definite lower bound on the radius for which this renormalization 
can be done. For low momenta, we find that the phase shifts are insensitive to the choice of regulator parameter, as 
one would expect since we have fit the effective range expansion up to second order. For radii close to the optimal 
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value, we observe that the agreement persists to higher momenta. Furthermore, for these nearly-optimal radii the 
quantity AT is small for momenta up to roughly 0.5 fin , thus indicating that the EFT expansion is under control 
for momenta below this scale. Again, we would not expect our EFT to be valid above this scale since it does not 
explicitly contain pions. 

We do find, somewhat disappointingly, that a straightforward effective range expansion of the data fits the actual 
phase shifts to considerably higher momenta than our optimal well. This indicates that the shape parameter for our 
optimal square well is considerably larger than the experimental shape parameter, which appears to be very small. 
While this situation could be remedied by adjusting our radius to obtain a better agreement with the data such an 
approach would not be systematic, as it would rely on the sensitivity of the phase shifts to the cutoff parameter R. 



III. AN EFFECTIVE FIELD THEORY WITH EXPLICIT PIONS: 

SQUARE WELL 



ONE-PION EXCHANGE AND A 



Naturally, the pionless EFT of Section || is only valid for momenta considerably less than m„. In order to construct 
a better EFT and so draw nearer to cases of interest in nuclear physics in this section we include pions in our EFT. 
At zeroth order the NN potential resulting from such an EFT Lagrangian is: 



Vfl(x) =C <4 3) (x) + V4(x), 



(3.1) 



where V n (x) is the 1 <So OPE potential, which, if we absorb the delta function piece into the contact interaction takes 
the form 



(3.2) 



where 
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167T/2 



(3.3) 



Strictly speaking at second order in the EFT we should calculate pionic corrections to this potential as well as 
including derivatives of the regulated contact interaction. However, as a first atte mpt we inc lude only these deriv ativ e 
terms, using as our second-order EFT interaction a sum of the interaction of Eqs. ( |2.3| ) and (2.4) and the OPEP fl3.2|). 

The radial Schrodinger equation in the I = partial wave is then: 



77 )^ 



1 

M 



2C 2 I -j- 

dr 



dU U{r) 



^ + ( \;ir) -! C^uir) 



d 2 5 R 
dr 2 



U{r) = EU{r). 



(3.4) 



The boundary co ndit ion at r = R is identical to the case with only a square well given in Eq. (2.6). 

The equation (3.4) cannot be solved analytically. Therefore, we must solve the differential equation numerically 
and impose a matching condition far from the origin in order to find the phase shifts. Once this is done, the analysis 
proceeds exactly as in Section JD. We note that because the pion potential is not compactly supportedjff is no longer 
clear that the results of Ref. [O] demand a lower bound on R. We will, nevertheless, see in Section |y] that such a 
lower bound does exist in this case. 

When we set C2 = and solve for the values of R and Cq which reproduce the 1 So NN scattering length and 
effective range, we obtain an optimal well width of 2.3 fm with Co = —3.3 fm 2 . A plot of the phase shifts as a function 
of k for this potential is shown in Figure f^. 



8 




k ( fm ) 



FIG. 7. A plot of phase shifts in radians versus momentum for the pion potential plus an optimal well. These are calculated 
from the numerical solution of Eq. (3.4) with R — 2.3 fm, Co = —3.3 fm 2 and C2 = 0. 



Next we renormalize the coefficients Co and C2 for other values of R. Figures || and || plot the renormalized values. 
We observe that both Co and C2 diverge near R = 1.4 fm, indicating that this is the smallest well that can be used 
to parameterize the NN interaction in the ^o channel when a one-pion exchange potential is included. 
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FIG. 8. A pl ot o f the renormalized values of the zeroth-order coefficient in the EFT expansion, Co, from the numerical 
solution of Eq. (3.4) that reproduces the low-energy 1 5'o NN scattering data, versus the radius, R, of the square well piece of 
the potential. 
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FIG. 9. A plot of the renormalized values of the second- order coefficient in the EFT expansion, C2, versus the radius, R, of 
the square well piece of the potential. 

Now that we have renormalized Co and C2, these values can be used to calculate the phase shifts for our second- 
order EFT potential. A comparison of the phase shifts for several values of R is shown in Figure [l0| We observe the 
same qualitative behavior as was seen in Figure ^. The phase shifts for all radii agree up to a point — as must be, 
since we have fit the first two terms in the effective range expansion. At this point, the results for potentials with 
radii far from optimal diverge rapidly. However, for values of R reasonably close to the optimal value, the agreement 
persists to a considerably higher momentum. 
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FIG. 10. A plot of the phase shifts in radians versus momentum for various radii from the numerical solution of Eq. 
The phase shifts for the optimal well are shown in bold. 



(3.4) 



Nevertheless, one sees that sensitivity to the particular regulator chosen creeps into the observables at a momentum 
of k k, 0.5 fm -1 (earlier if the well is ridiculously large). This is not surprising since this is roughly the momentum 
at which details of the well structure begin to be probed. One might think that this is an argument for removing the 
regulator from the problem by taking R to zero. However, in this work we have shown that R cannot be decreased 
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below about 1.4 fm. Thus, we appear to be forced to use a regulator which limits the applicability of the EFT to 
momenta below some maximum, A. This maximum lies considerably below the expected range of validity of the EFT. 
This is not a desirable situation. Similar comments apply to the results of Fig. ||, but there the pion is integrated 
out, and so A is of roughly the same magnitude as the expected maximum momentum for the EFT. What we see in 
Fig. |h] is that including the pion explicitly does not greatly decrease the sensitivity of the results to the radius of 
the short-range potential. This occurs because the large effective range of this problem means that the "short-range" 
potential cannot really be made short range at all, but instead must extend out to at least 1.4 fm. This suggests the 
possibility of a breakdown in the scale separation which is essential for power-counting arguments to be applicable. 

Figure |ll| compares the phase shifts for the OPE potential plus optimal well with the Nijmegen np partial wave 
analysis for the 1 S'o channel pj^ . This plot shows good agreement up to the pion production threshold. The divergence 
at higher momenta is sharp but not unexpected since calculated phase shifts are positive for all momenta while the 
experimental phase shift becomes negative just above pion threshold. This behavior could be reproduced if we tuned 
the potential specifically to do so. 
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FIG. 11. The solid line in this plot is the phase shift as a function of momentum for the OPE potential and an optimal 
square well. It is compared with data points from the Nijmegen np partial wave analysis in the 1 So channel M] 



However, we observe that there is really no reason to compare the Nijmegen phase shifts with the results from the 
optimal well. The phase shifts from any of the other wells whose phase shifts are plotted in Fig. |l^ can be compared 
to the data too. Upon doing this it is clear that the optimal well happens to be a potential whose structure affects 
the phase shifts in a way which brings them closer to to the experimental data. In other words, the good agreement 
seen in Fig. ^ is fortuitous, and is dependent upon details of the particular regulator used. We cannot expect such 
good agreement in general. 

Finally, we once again compare the the zeroth- and second-order t-matrices for different values of R. We will 
examine the fractional difference ATn(k) as defined in Eq. ( 2.14 ). The second-order t-matrix, (k), is calculated 



using the pion potential and a square well with second derivative terms. As in Section [n| the zeroth-order t-matrix is 
calculated with only the zeroth-order potential ( |3.l| ) , fitting the scattering length to match that of the NN interaction 
in the x So channel. The results are plotted in Figure [l2| Again the behavior is very similar to the case where the 
pion was integrated out. The inclusion of the pion has not decreased the size of this quantity significantly, because 
the "short-range" square well is itself of roughly one-pion range. 
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FIG. 12. A plot of the fractional difference between the zeroth- and second-order t-matrices versus momentum for several 
values of R. For an optimal well, with R—2.3 fm, this quantity is identically zero. 



IV. A LOWER-BOUND ON THE RANGE OF THE SHORT-RANGE INTERACTION 



In the previous section we saw that the radius of the square well potential could not be reduced below 1.4 fm if the 
scattering length and effective range generated by the sum of a short-range and one-pion exchange potential were to 
agree with the experimental results. In fact, there exists a lower bound on this radius which is independent of the 
choice of regulating potential. In this section we extend the arguments of Ref. |l2| in order to show that, if the theory 
is to fit the scattering length and effective range, then there is a general absolute lower bound on the range of the 
non-one-pion-exchange piece of the interaction. 

Consider the radial Schrodinger equation for S-wave scattering in the case where both a one-pion exchange potential 
and an arbitrary (possibly non-local) short-range potential, Vr are present: 



1 d 2 u E (r) 
M dr 2 



dr'rVff(r, r')r' ue{t') 



-u E {r) = Eu E {r), 



(4.1) 



with ue(0) = 0. Here Vr{t : r') = for r > R or r' > R. This, of course, was the case discussed in the previous 
section for certain specific Vrs. However, in this section, instead of enquiring as to the exact nature of Vr, we now 
consider a solution, ve (r) of the equation 



1 d 2 v E (r) 
M dr 2 



-v E (r) = Ev E (r). 



(4.2) 



ve is chosen so as to match onto the asymptotic wave function, u E {r) — sin(/cr + 5 E ), with k = v 'ME and Se the 
expe rimental phase shift, as r —* oo, and is normalized so that «e(0) = 1. Suppose now that ue(t) is a solution of 
(4.1), with the parameters of Vr adjusted so that u E has the experimentally observed asymptotic behavior. Suppose 
also that ue is normalized so that it agrees with ve{t) at r = R. Given this normalization the two wave functions 
agree on [R, oo). However, they differ on [0, R], in that u E (0) = 0, while Ue(0) = 1. 
Now, going through the arguments displayed in Jl^| yields: 



dv 2 
dr 



r=0 



dv\ 
dr 



= {k 2 



r=0 



k 2 ) I dr [v2(r)v±(r) - u 2 (r)ui(r)], 



(4.3) 



where v 2 and v\ are obtained by solving (4.2), with the a ppro priate boundary conditions at r — R, at two different 
energies E 2 and E%. Similarly, u 2 and «i are solutions to (|4.1|) at the same two energies. From Eq. (4.3) we find 



d ( dvE 
dE \~dr~ 



r=0 



M 



dr [v E {r) - u E {r)\ 



(4.4) 
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Since UE{r) may be chosen to be real, and ue and ve agree for r > R, it follows that 



d ( <1ve 

dE \~dT 



< M 



dr v%(r) 



(4.5) 



Once the asymptotic behavior of ve is specified the function ve(t) is independent of R, depending only on the 
experimental phase shift Se- By contrast, the right-hand side is a function of i?, but only through the integral's upper 

I ) is positive, then as R is decreased towards zero a value of R will be reached for which 



bound. Thus. 



if 

11 dE 



(4.5) will be violated. 

In the case where there is no one-pion exchange interaction the wavefunction ve{t) is 



ve{t) 



sin(fcr + 5 1 
sm(8 E ) 



(4.6) 



and so, Eq. (|4.5|) becomes: 



dE 



(fccotfe) < 0, 



(4.7) 



as claimed in the Introduction. As mentioned in Section ||, this may be derived from an old result of Wigner Jl5| . 

In the case where one-pion exchange is included the function ve{t) may be calculated numeri cally given experimental 
phase shift data. When this is done using the experimental values for a and r e given in Eq. ( 2.12| ) the result: 



_d_ 

dE 



dvE 
dr 



= 10.0, 



(4.8) 



r=0/ E=0 



accurate to two significant figures, is obtained. (In practice this is best done by fixing the logarithmic derivative at 
some small, but finite, distance, integrating out to large distances and m atch ing to the experimental phase shifts, and 
only then integrating in to the origin.) We may now check whether Eq. ([h^) is violated at E = 0. In Fig. [l^ we plot 
the right-hand side of Eq. (p~5|), taken at E — 0, 



f(R) = M drvl{r), 



(4.9) 



as a function of R. The constraint ( [4.5[ ) is violated once R < 1.1 fm. It follows that any short-range potential of range 
less than 1.1 fm which is used in the Schrodinger equation (4.1) will not be able to fit the experimental scattering 
length and effective range. 
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FIG. 13. A comparison of the two sides of Eq. ( k.5L e valuated at E = for NN scattering in the 1 5o channel. Observe that 
the function f(R) is less than the numerical bound (|l.q) for radii less than 1.1 fm. 
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V. CONCLUSIONS AND DISCUSSION 



In this paper we have shown how a finite-range potential can parameterize the short-range physics in an effective 
field theory (EFT) approach to the nucleon-nucleon interaction. We have chosen a square-well potential since it allows 
a simple and clear analysis using elementary quantum mechanics. Other forms of potentials could be used and should 
have little effect on the results. We have considered both the case where all exchanged particles are integrated out, 
and the case where a one-pion exchange potential is retained. 

We have done EFT calculations at both zeroth and second order. In doing so we have adopted the approach of 
Weinberg, and done the power counting in our EFT in the potential, rather than in the amplitude. The zeroth-order 
calculation then involves a potential which is the sum of a square well and one-pion exchange (if explicit pions are 
present), while our second-order EFT calculation adds derivatives of the regulated contact interaction. (It should 
be noted that our "second-order" EFT calculation in the theory with pions ignores some two-pion exchange graphs 
which are, strictly speaking, of the same EFT order as graphs we have included in our potential. However, we believe 
the inclusion of these graphs in the calculation will not qualitatively alter our conclusions.) For a given regulator 
parameter the coefficients in the short-range interaction are fitted to the 1 5'o scattering length and effective range. 

An optimal choice for the well radius was found in both the theory with explicit pions and that without. At that 
well radius the second-order piece of the potential is identically zero and yet the 1 5'o scattering length and effective 
range are reproduced. When one attempts to renormalize the coefficients in the EFT expansion so as to fit the 
experimental data at other well radii, there turns out to be a lower bound on the well radii for which this can be done 
successfully. This suggests that this method of regulating the theory does not allow one to reach the limit of truly 
zero-range interactions. However, it should be noted that other, non-equivalent, regulation schemes may allow one to 
define a Schrodinger equation containing contact interactions. 

For potentials with radii close to the optimal value the second-order corrections to the t-matrix are small, and the 
phase shifts are similar to those produced by the optimal well. This lack of sensitivity to the regulator parameter 
indicates that it might be possible to develop a systematic power-counting scheme along the lines of the calculations 
performed in this paper. However, our results show that such a power-counting scheme can only hope to be successful 
for momenta up to about 100-150 MeV. This is a range of validity which excludes many interesting nuclear physics 
applications. It is somewhat surprising that the cutoff cannot be reduced beyond 1.4 fm for the case of a square- well 
regulator. Indeed, the results of Section [V show that, no matter what regulator is used, it must extend beyond 
r = 1.1 fm. This means that the physics which, together with one-pion exchange, explains the effective range in the 
1 S'o channel, is not particularly short-ranged. 

In the case where the OPE potential was included, we found that the optimal well fits the experimental phase shifts 
well. However, if other regulator parameters are chosen the data is not fit particularly successfully beyond about 
0.5 fm '. 

There are a number of points to observe in comparing the approach of this paper with that recently advocated by 
Kaplan et al. Q. First, we note that the regulation of the delta function interactions used in that work was completely 
different to that employed here. Whereas we have defined a contact interaction as the limit of a sequence of square 
wells of decreasing radii, Kaplan et al. have calculated sets of loop diagrams using true contact interactions and then 
renormalized the resulting infinities using dimensional regularization and the MS scheme. The two approaches to 
regulating the delta function potential are not equivalent. Specifically, in our approach the range of the short-range 
potential cannot be taken to zero if the renormalization conditions on its coefficients are to be satisfied. 

Second, upon keeping the short-range potential of finite range, we see that using the Weinberg approach to power 
counting allows the explanation of both the scattering length and the effective range in an EFT which is valid to 
momenta k « 0.5 fm -1 . This is in contrast to the results of Kaplan et al. who found that, given their method of 
regulation, if the scattering length and effective range were to be explained in a Weinberg power-counting approach 
then the domain of validity of the resulting EFT was very small. It was this that led Kaplan et al. to define a 
new power-counting scheme in which power counting was employed for the inverse scattering amplitude. With the 
form of regulation used in this paper power counting may always be applied to the potential without a poor radius 
of convergence for the EFT resulting. Of course, our effective field theory expansion is open to question since our 
regulator was never, and indeed can never be, removed from the theory. However, we have demonstrated that when 
the regulator parameter is kept finite and within sensible bounds, there is not great sensitivity to it for processes 
involving momenta up to about 100 MeV. A formal power-counting scheme for an EFT approach where the range of 
the short-range interaction is always kept finite remains to be worked out. 

Finally, from the standpoint of correctly reproducing the experimental data, in the case without pion exchange the 
approach of Ref. j|] (which there is identical to effective range theory) does considerably better than our "optimal 
well". When the pion is included, the two approaches appear to yield nearly identical results, although the lack of an 
a priori reason for choosing the optimal well over wells of other radii means that we cannot unambiguously say that 
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the phase shifts predicted by our approach are those shown in Fig. ED]. Above 100 MeV the predictions from potentials 
in which different regulator parameters were used differ considerably. This is a troubling result, and casts doubt on 
the efficacy of the approach discussed here. However, strictly speaking, the momenta to which the EFT of Kaplan et 
al. applies is of roughly this size, although a comparison of their results with experimental data indicates that good 
agreement is obtained to much higher momentum than 100 MeV. Thus, it remains to be seen if either approach can 
be successfully used in describing other nuclear physics problems. 
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